Quasi-equilibrium models for triaxially deformed rotating compact stars 
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Quasi-equilibrium models of rapidly rotating triaxially deformed stars are computed in general 
relativistic gravity, assuming a conformally flat spatial geometry (Isenberg- Wilson-Mathews for- 
mulation) and a polytropic equation of state. Highly deformed solutions are calculated on the 
initial slice covered by spherical coordinate grids, centered at the source, in all angular directions 
up to a large truncation radius. Constant rest mass sequences are calculated from nearly axisym- 
metric to maximally deformed triaxial configurations. Selected parameters are to model (proto-) 
neutron stars; the compactness is M/R = 0.001,0.1,0.14,0.2 for polytropic index n = 0.3 and 
M/R = 0.001,0.1,0.12,0.14 for n — 0.5. We confirmed that the triaxial solutions exist for these 
parameters as in the case of Newtonian polytropes. However, it is also found that the triaxial 
sequences become shorter for higher compactness, and those may disappear at a certain large com- 
pactness for the n = 0.5 case. In the scenario of the contraction of proto-neutron stars being subject 
to strong viscosity and rapid cooling, it is plausible that, once the viscosity driven secular instabil- 
ity sets in during the contraction, the proto-neutron stars are always maximally deformed triaxial 
configurations, as long as the compactness and the equation of state parameters allow such triaxial 
sequences. Detection of gravitational waves from such sources may be used as another probe for the 
nuclear equation of state. 



I. INTRODUCTION 

Rapidly rotating compact objects are expected to 
be formed as new born neutron stars after stellar 
core collapses, or as differentially rotating hypermas- 
sive neutron stars after binary neutron star merg- 
ers. Accretion onto neutron stars in X-ray binaries 
can also lead to rapid rotation. All of these have 
been extensively studied as strong sources of gravita- 
tional waves for the ground based laser interferometers 
LIGO/GEO600/VIRGO/TAMA (See e.g. [1 and refer- 
ences therein). 

Classical models of rotating stars are a class of ellip- 
soidal figures of equilibrium; self-gravitating rotating in- 
compressible fluids in Newtonian gravity. Such solutions 
include sequences of axisymmetric Maclaurin ellipsoids, 
or non-axisymmetic Jacobi, Dedekind, or Riemann S- 
type ellipsoids [2I]. These models are used to study the 
secular evolutions of rapidly rotating stars due to the 
viscosity and the radiation back-reaction of gravitational 
waves 01 • Lai and Shapiro [1| have developed an ellip- 
soidal approximation to the rotating polytropes, and ap- 
plied the model to clarify the secular evolution of rapidly 
rotating neutron stars in detail, and more recently fo- 
cused on the viscosity driven secular instability Q. 

As discussed in P, d, and shown by a number of 
numerical simulations of rapidly rotating compact stars, 
core collapse, and binary neutron star mergers @, 0, @, 
a transient triaxially deformed compact object may sur- 
vive within a secular time scale. In this paper, we con- 
sider uniformly rotating models of such triaxially de- 
formed compact objects, an extension of the Jacobi el- 
lipsoid in general relativity. In Newtonian gravity, such 
solutions exist for rotating polytropes with polytropic in- 
dex n < 0.808 [9i]. Here, the polytropic equation of state 



(EOS) p = Kp^~^^^" relates the pressure p with the baryon 
rest mass density p. 

In general relativity, such configurations are not in 
equilibrium due to the back-reaction of gravitational ra- 
diation. However as in the case of quasi-equilibrium ini- 
tial data of binary neutron stars, triaxially deformed uni- 
formly rotating stars are in quasi-equilibrium, as long as 
the gravitational luminosity is small enough that the en- 
ergy radiated away within a rotational period is small 
compared to the binding energy of the star, which is al- 
ways the case, and if the viscosity is strong enough for 
the flow field to become uniformly rotating during the 
evolution. Therefore, as an important application, a se- 
quence of uniformly rotating quasi-equilibrium solutions 
may model a secular evolution from the proto-neutron 
star to the neutron star in the strong viscosity limit, and 
each solution may serve as the initial data for the general 
relativistic hydrodynamic simulations of such objects. 

Models of rapidly rotating neutron stars have been 
extensively studied as stationary, axisymmetric, perfect- 
fluid spacetimes while less attention has been paid 
to uniformly rotating triaxial solutions, not only because 
they are not exact equilibria due to gravitational radia- 
tion reaction, but also because, in early calculations [ll[, 
models of the EOS above nuclear density did not allow 
large enough values of r/|iy| ~ 0.14 where a triaxial 
sequence is expected to bifurcate from an axisymmet- 
ric sequence. Here, is the ratio of kinetic en- 
ergy T to gravitational potential energy W. However, 
as seen, for example, in [1^], a value n ^ 0.5 (an effec- 
tive adiabatic index F ~ 3) may be possible for recent 
models of EOS for high density nuclear matter, above 
P > Pnuc ~ 2 X 10^'^g/cm'^, and, as we will see below, tri- 
axial quasi-equilibrium solutions do exist even in strong 
gravity for relatively small polytropic indexes, such as 
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n = 0.5 or n = 0.3, as in the Newtonian case. 

About a decade ago, Nozawa succeeded in comput- 
ing uniformly rotating triaxial polytropes in general rel- 
ativistic gravity in his thesis flS], although his cal- 
culations were limited by the computational resources 
to low resolutions. A few studies approximating the 
fluid as an ellipsoidal configuration in general relativistic 
gravity have been made and perturbative analyses 
have located the bifurcation point suggesting the exis- 
tence of solutions having triaxial bar-mode deformations 
[H [H [13, E, [llll^l ■ Our computations of triaxially de- 
formed stars can be used to locate the instability points 
on the axisymmetric sequence as well as to estimate the 
gravitational wave amplitude and luminosity from such 
objects. 

In this paper, we present our first results on triax- 
ial configurations of rapidly rotating general relativistic 
stars as models of neutron stars in extreme rotation. We 
assume a conformally flat spatial slice, and solve the 
constraints and spatial trace of the Einstein equation 
(Isenberg- Wilson-Mathews (IWM) formulation) Hi] 
^ . This is different from the formulation used by Nozawa 
[T^ l in which the line element is chosen to be the same 
form as that of stationary axisymmetric spacetime, but 
an azimuthal dependence is allowed (see [2^ for the same 
formulation). The formulations and the code are de- 
scribed in the next section. For testing the code, selected 
axisymmetric solutions are compared with the results in 
the literature, and the bifurcation points of axisymmet- 
ric and triaxial sequences in weak gravity are examined. 
Then, we present the results of deformation sequences 
of constant rest mass systematically in the range of two 
parameters, the compactness M/R and the polytropic 
index n, appropriate for realistic neutron stars. Appli- 
cations of such triaxial solutions in the contraction of a 
newly born proto-neutron star are briefly discussed in 
the final section. Throughout the paper, we use units 
such that G = c = 1. For our tensor notation, we adopt 
the use of Greek letters for spacetime indices, and Latin 
letters for spatial indices. 



to the timelike vector t", which is tangent to a curve 
t {t,x), X € E, by = an°'+/3'^, where a is the lapse, 
and f3" the shift which satisfies (3°'na = 0. Restricting 
the projection tensor jap = gap + naUp, a spatial metric 
7a6(t) is defined on Ej. In the IWM formulation, the spa- 
tial metric is chosen to be conformally flat, 'jab = ip'^fab, 
where fat is a flat metric on each slice, and is a con- 
formal factor. Then the metric gap takes the form 

ds"^ = -a^dt^ + tl/fij{dx' + /3'dt)(dx^ + (3^dt), (1) 

in a chart {t, a;*}. The extrinsic curvature of the foliation 
is defined by 



Kab ■= --Z-dtJab + -Z-i^Plab- 

2a 2a 



(2) 



In the IWM formulation, the Einstein equation is de- 
composed with respect to the foliation, and the following 
5 components Eqs.®-® are chosen to be solved for the 
five metric coefficients {"f/;, a, on the initial slice Y^q: 

(Gap - SnTap) n'^nP = 0, (3) 
{Gap - SnTap) ja^n^ = 0, (4) 

{Gap - SnTap) (7"'' + = 0, (5) 

where the first and second equations are the constraints. 
These equations are written in the form of elliptic equa- 
tions with the non- linear source terms, respectively, 

= (^AabA-' + ^K^^ - 2^7AVh, (6) 

o~ loo- o 4 o 

A(3a + -DaDbp^ = -2aAa''Db\n^ + -aDaK 
6 a 6 

+ l&TTaja, (7) 

+a'<P^ {^-AabA'^'' + ^K^\ + 2Tia'4,HpY^ + 25), (8) 



II. FORMULATION AND NUMERICAL 
METHOD 

A. IWM formulation 

The IWM formulation for computing spatially con- 
formally flat initial data we use for computing non- 
axisymmetric quasi-equilibrium of rotating compact star 
is briefly described. The spacetime M. = Rx E is foliated 
by the family of spacelike hypersurfaces Et — {<} x E. 
The future-pointing timclikc normal n" to Ej is related 



^ The vaUdity of the IWM formulation for axisymmetric configu- 
rations is discussed in l23l 12411 



where K is the trace of Kab, Aab its tracefree part, Pa 
the conformally weighted shift deflned by fia = fabP^ and 

~ o o 

pa _ ^ A is the flat Laplacian and Da is the covariant 
derivative with respect to the flat three-metric fab- The 
source terms of matter are defined by '■— Tapn^n^, 
ja ■— —Tapja^n^ , and 5* := Tapj"^ ■ Expressions of 
the sources in terms of the metric potentials and fluid 
variables are given in Appendix [XI 

We choose a maximally embedded slice K = — dtK. 
Because the spatial metric is conformally flat, Aab does 
not involve time derivatives of the spatial metric. 



Aab — ^ (^-pfab — -^fabf'^'^'i^pfa 



(9) 



lb* / o -. o -. 2 o\ 
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where £^ denotes the Lie derivative with respect to (3°". 
The field equations Eqs.®-® are thus rewritten 



(11) 



a 
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A{a^) = -a^^AabA'^' + 2^a^^{p^ + 2S). (13) 

o 

Eq. is decomposed further to improve the accuracy 
in numerical computation. Following the decomposition 
proposed by Shibata [2^ , we write Eq. as 



and introduce 



(14) 



(15) 



o 

where a;" are coordinates that satisfy DaX = Sa . Sub- 
stituting the decomposition (fT5|) into Eq. fTl)) yields 

A;3, + ii3aA/3'' = ABa + ^Da{AB - J" ABb) = Sa. 

(16) 

O 

The elliptic equations ABa = Sa and AB — x ABb = 
are separated, and the former is substituted to the latter: 

ABa = Sa := -2aA,^i3bln— + 16™ia, (17) 

a 



AB = x'^Sa. 



(18) 



The potentials {i?a, -B} are solved for simultaneously, and 
the shift Pa is reconstructed from Eq.ljTS]). 



B. Formulation for a relativistic fluid in 
equilibrium 

A perfect fluid is described by the stress-energy tensor 



Tap = (e + p)UaUp + pgap, 



(19) 



where is the 4- velocity of the fluid, p its pressure, and 
e the energy density. As a consequence of the Bianchi 
identity, the stress-energy tensor is covariantly conserved: 



VpTj = 



(20) 



When the fluid is close to equilibrium, one can obtain 
a simpler set of equations. Introducing the specific en- 
thalpy defined by h := (e + p)/p, where p is the baryon 
rest mass density, the left hand side of Eq. (PD|) can be 
written 



/3 _ 



+ huaVp{pu'^) - pTVaS, (21) 



where s is the specific entropy. In the derivation, the 
local first law of thermodynamics dh = Tds + dp/p was 
used. In local thermodynamic equilibrium, one can also 
assume the conservation of baryon mass. 



Vo(pw") 
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:£«(P\/^) = 0. 



(22) 



Consequently, the conservation of specific entropy along 
the fluid world line, 

u'^VaS = £uS = 0, (23) 

and, the relativistic Euler equations, 

u^^V p{hua) + V ah = £u{hua) + Vah = 0, (24) 

are obtained. Assuming the flow field to be isentropic 
everywhere inside the neutron star matter, s = const, we 
have a one-parameter equations of state (EOS) p — p{p)- 
We assume a stationary state in the rotating frame for 
the fluid source. Imposing symmetry along the helical 
vector k" — t" + $10" where f2 is a constant angular 
velocity of a rotating frame, we have 



£k{pu*V-9) = 0, and ^a£k{hua) = 0, 



(25) 



with u* interpreted as the scalar u'^V at- For a corota- 
tional flow, u" = M*fc", the rest mass conservation be- 
comes trivial, and the relativistic Euler equation has the 
first integral 



£ 



constant, 



(26) 



where £ is the injection energy. From the normalization 
of the four velocity UqM™ = — 1, one obtains 



1 



(27) 



where w° = -I- ^cf)"- . 

As a first step in the calculation of a highly deformed 
triaxial compact star, we assume a simple polytropic 
EOS, 



p = up 



1+1/ti 



(28) 



where k is a constant, and n is the polytropic index. 
Then h is related to p/p by 



l + (n+l) 



P 
P' 



(29) 



We also refer to the polytropic exponent F defined by 
F := l-f 1/n. 



C. Numerical computation 

The Poisson solver and the iteration scheme used to 
solve the system of elliptic equations with non-linear 



4 



source terms, are similar to the ones used in a previ- 
ously developed initial data code for binary black holes 
and neutron stars (27l. [28l|. However, the code itself has 
been completely rewritten, so that further extensions can 
be incorporated easily. One of the revisions of the code 
is that no symmetry is a priori assumed on the spatial 
slice Eq; that is, the spherical coordinate grids centered 
at the source cover all angular directions, up to a certain 
large truncation radius. Hence, for example, asymmet- 
ric magnetic fields may be later included without major 
modifications to the code. Computation of binary solu- 
tions using the same coordinate setup is also possible. 
The other major change is a simpler, more robust choice 
of finite differencing. In this section, we briefly describe 
the necessary steps for constructing the code, which are 
1. Spherical coordinates and the length scale, 2. Sum- 
mary of variables and equations for coding, 3. Poisson 
solver, 4. Grid spacing, 5. Finite differencing and itera- 
tion, 6. Computation of a sequence of solutions. 

1. Spherical coordinates and the length scale 

The slice So is covered by a spherical coordinate patch 
(r, 6*, (j)) g [ra, Tb] X [0, tt] x [0, 27r]. For a single star calcula- 
tion, the radial coordinate extends from the center of the 
star r = ra = to the asymptotic radius r = = 10''i?0i 
where Rq is the radius of the neutron star along the semi- 
major axis, defined by the 9 = tt/2 and (/) = 0, tt lines. We 
also refer to Cartesian coordinates {x, y, z) whose pos- 
itive X, y and z directions are along (0,0) — (7r/2,0), 
(7r/2, 7r/2) and = 0, respectively. 

The quantity i?o is introduced as an additional param- 
eter in the formulation used in our code, normalizing the 
radial coordinate as 

f = r/i?o. (30) 

For a polytropic EOS, one can rescale the length using 
the polytropic constant 

i?Oi or simply setting 
K= \ (see e.g. [l^)- As a result, we have three parame- 
ters {ri,£,i?o} in our formulation. 

Furthermore, we introduce surface fitted coordinates 
on which the fluid variables are deflned. Assuming that 
the surface of a neutron star can be described by a func- 
tion of the angular coordinate R{9,(j)) = RoR(9,(t)), the 
surface fitted coordinates (ff , 9f, (j)f) are defined by 

ff :=f/i?(0,0), 9f = 9, 4>i^(t>, (31) 

where ff is deflned in a region ft G [0, 1]. 

2. Summary of variables and equations for coding 

As mentioned in Sec. Ill Al the fleld equations (fTTj) . 
([13]), P7)) and are solved for the metric potentials 
{■0, aV', . B} and, as in Sec. IH Bl and lll C 11 a comoving 
fluid in equilibrium is characterized by one fluid variable. 



which is chosen to be the relativistic enthalpy {/i}, and 
three parameters {n, f , i?o}- 

The fleld equations are normalized to have the fol- 
lowing form; representing each of the metric potentials 
{V',aV',S",^} by 

A$ = 5g + (32) 

o 

where the fiat Laplacian A corresponds now to the nor- 
malized coordinate f . The source term includes the 
metric potentials and their derivatives, while S-ai also in- 
cludes the matter variables and the parameters 
while the dependence on the length scale i?o is explicitly 
separated in Eg. ([5^ . 

The fluid variable {h] is determined by Eg. ipSl) cou- 
pled to the EOS dSHD, and the relations ^ and (1^ . 
The three parameters {r2,f,i?o} are determined by the 
following three quantities: the surface radii along two of 
the three semi-major axes, and the value of the central 
density. These quantities are used to impose three con- 
ditions on Eg. ([26]) . which are solved with respect to the 
three parameters {i7,£,i?o} in each iteration cycle. 



3. Poisson solver 

The elliptic equations (HH), (UHl) (HZ]), and ^ are inte- 
grated on the spherical grid using Green's formula. Rep- 
resenting each of the potentials {^, atj}, -B°, B} by 3>, the 
latter is given by 

^{x) = -— / G{x,x')S{x')Sx' 

47r Jy 

+ 1 j [G(:r,:r')V'"ci>(x') 

-^{x')V"'G{x,x')]dS'^. (33) 

where x and x' are positions, e y C Sq. We choose 
the Green function G{x,x') without boundary, 

AG{x,x') ^ -AttS{x-x'), (34) 

and perform a multipole expansion in associated Legen- 
dre functions, 

G[x, X = = yge[r,r) > e™ -— -r 

X " (cos 9) Pf " (cos 9') cos m(v9 - ^'), (35) 
where the radial Green function gi{r,r') is defined by 

9i(.r,r') = r> := sup{r,r'}, r< := inf{ r,r'}, 

(36) 

and the coefficients e™ are equal to ep = 1 for m = 0, 
and e™ = 2 for m > 1. 
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Ta : Radial coordinate where the grid Ti starts. 

ri, : Radial coordinate where the grid ends. 

rc : Radial coordinate between ra and rt, where the 

grid changes from equidistant to non-equidistant. 

Nr : Total number of intervals Ar; between Ta and r^. 

rir : Number of intervals Ar^ between ra and rc. 

Nf : Total number of intervals Af^ for f G [fa, R{0, (j))]. 

Ng : Total number of intervals A9i for 9 £ [0,7r]. 

A^,^ : Total number of intervals A<j>i for (p € [0, 27r]. 

TABLE I: Summary of grid parameters. 

4- Grid spacing 

The field equations in the integral form ([55]) are dis- 
cretized on the spherical grids, and iterated until con- 
vergence is achieved. Our code allows us to use any 
non-equidistant grid spacing in all the spatial coordi- 
nates, {ri,9j,4)k), i = 0, •••,7Vr, j = 0, •••,7Ve, and 
k — 0, • • • , N^. The radial grid points are equidistant 
in the region [ra,rc] and non-equidistant in [rc,rf,], as 
follows: writing Ar,; :— Vi — r.i_i, we have 

Ar, = Ar = Is^ll^ for i = 1, . . . , n^, (37) 
Ari = fcAr.,_i for i = n,. + 1, • • • , Nr, (38) 
where the constant k is determined from the relation 

n-rc = Ar . (39) 

k — \ 

We choose equidistant grid spacing for 9j and 0^, that 
is, AOi = AO = ir/Ng, and A(j)i = A0 = 2it/N^. Our 
notations for the grid points are summarized in Table U 

5. Finite differencing and iteration 

For the numerical integration of Eq. we select the 
mid-point rule. Accordingly, source terms are evaluated 
at the middle of successive grid points. The linear inter- 
polation formula and the second order Lagrange formula 
are applied for computing the source term fields and their 
derivatives respectively, at the mid-points of the r, 6 and 
4> grids. 

The reason for selecting a rather low (second) order 
finite difference scheme is the following: When the field 
quantities vary rapidly, such as at a density discontinuity 
in a neutron star, higher order interpolating formulas as 
well as finite difference formulas tend to overshoot near 
the region, and may cause a non-convergent iteration. 
To overcome this behavior, one may either (i) separate 
the computing regions at the discontinuity, or (ii) use 
lower order polynomial approximations. With the first 
approach, pseudo-spectral methods have been success- 
fully implemented by [2!| and achieved an evanescent 



error. We select the second idea to keep the code as sim- 
ple and flexible as possible, and improve the accuracy by 
simply increasing the number of grid points. 

In each iteration cycle, the Poisson solver (|33p is called 
for each variable. Writing the L.H.S. of Eg. ([55)1 as $, 
each field variable is updated from the iVth iteration cycle 
to the {N -t- l)th in the manner 

$(^+1) = aI> + (1-A)$(^\ (40) 

where the softening parameter A is chosen to be 0.3 ~ 0.5 
for accelerated convergence. Then we check the relative 
difference of successive cycles 

±Ll 1. L (41) 

I + \ ^{N) I' ^ > 

as a criteria for the convergence. We typically stop the 
iteration when this quantity becomes less than 10~^. 

The method used in this code may be considered as an 
extension of the one developed by Ostriker and Marck 
(1968) [13] for Newtonian rotating stars, and by Ko- 
matsu, Eriguchi, and Hachisu (1989) for relativistic ro- 
tatin g (a xisymmetric) neutron stars, known as the KEH 
code [Sli. 

6. Constructing a sequence of solutions 

We compute constant rest mass sequences of isentropic 
equilibrium solutions in the IWM formulation. Constant 
entropy is modeled by setting the parameter k, in the EOS 
to a constant. In a quasi-equilibrium evolution of a ro- 
tating star, the angular velocity remains constant when 
the viscosity of matter is dominant. Then the solution 
sequence approximately models an evolution, with an er- 
ror that includes neglecting the increase of entropy due 
to viscosity. 

Each solution is computed by setting the central value 
oi q :— p/p to q — qc, the ratio i?(0, <j>)/Ro of the stellar 
radii along the z and x axes for the axisymmetric con- 
figuration, and the ratio i?(7r/2, 7r/2)/_Ro of the y and x 
axes for the triaxial configuration. To compute a con- 
stant rest mass sequence, one iterates Mo{qc) changing 
qc until Mq converges to the specified value. We use a 
discrete Newton-Raphson iteration for the rest mass. 

A sequence of rotating star solutions with a certain 
EOS is labeled by the compactness M/R of a non- 
rotating spherical star having the same rest mass Mq. 
We denote the gravitational mass of this spherical star 
by M, and the compactness by M/R. This labeling for 
each sequence is possible as long as it is a normal se- 
quence that has a stable spherical star in the limit that 
the angular velocity Q goes to zero, which is not the case 
for a supermassive sequence. In the following, we focus 
on the normal sequences whose compactness is close to 
its value for a neutron star, around M/R ~ 0.1 — 0.2. 

Formulas used for computing the rest mass Mq, ADM 
mass Madm, Komar mass Mk, as well as the total an- 
gular momentum J are presented in Appendix [X] For 
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Type rg rb Vc Nr n,- Nf- Ne L 

I-l W 1.25 60 20 16 24 48 12 

1-2 10'' 1.25 90 30 24 36 72 12 

1-3 10" 1.25 120 40 32 48 96 12 

1-4 W 1.25 180 60 48 72 144 12 

I- 5 10'' 1.25 240 80 64 96 192 12 

II- l 10* 1.25 120 40 32 24 48 8 
II-2 10* 1.25 180 60 48 36 72 10 
II-3 10* 1.25 240 80 64 48 96 12 

TABLE 11: Coordinate parameters, and the number of grid 
points with diflerent resolutions. L is the highest multipole 
included in the Legendre expansion. 



polytropic EOS, one can normalize these quantities by 
a certain power of the polytropic constant shown 
in the same Appendix. Hence we choose k — I units to 
present solution sequences. 

A sequence of solutions with constant rest mass is con- 
sidered as an evolutionary track of adiabatic changes in 
quasi-equilibrium. Under this assumption, the solutions 
in each sequence are parameterized by the angular veloc- 
ity fl, and the first-law relation 



e ft Madm T/\W\ 6c 

Present 0.4614 0.5252 0.1247 0.04281 0.7911 

CST 0.4592 0.5232 0.1247 0.04253 0.7911 

Present 0.6370 0.6672 0.1264 0.08711 0.6613 

CST 0.6360 0.6658 0.1266 0.08705 0.6613 

Present 0.7581 0.7222 0.1281 0.1310 0.5614 

CST 0.7585 0.7214 0.1284 0.1314 0.5614 



TABLE III: The numerically obtained values are compared 
with those based on Table 1 of CST. Model parameters of 
these solutions are n = 0.5, M/R = 0.298, Mo = 1.484 x 10"\ 
and M = 1.230 x 10"\ The quantities from CST are inter- 
polated to have the same central energy density using the 
four-point Lagrange interpolating polynomials. The gravita- 
tional potential energy W is defined by Eq. (IA9p . The eccen- 
tricity e is defined by e := \/l — {Rz/RxY where the radii 
Rx , Rz along the x and z axes are measured in proper length 
as in Eq. (|A11|I . 

number of multipoles L as shown in Table HIl ^. 

When a sufficient number of multipoles is kept, the dif- 
ferences between numerically computed quantities with 
different resolutions and their exact value are written 



SMabm = ^SJ (42) 
is satisfied, as proved in [s^ . 

III. CODE TEST 
A. Axisymmetric solutions 

Axisymmetric solutions calculated by our new code 
are compared with the results in the literature (23l . [s^ . 
We show the results of comparisons for models presented 
in Table I of Cook, Shapiro and Teukolsky [23] (here- 
after CST), which correspond to a solution sequence with 
constant rest mass Mq — 0.14840 for the case with the 
polytropic index n = 0.5. This value of Mq is close to 
the maximum rest mass of a non-rotating spherical so- 
lution; the gravitational mass and the compactness of 
the same non-rotating solution are M = 0.12304 and 
M/R = 0.29605. In Table [ml selected solutions calcu- 
lated with the highest resolution 1-5 in Table [Til are com- 
pared with the results shown in Table I of 23] . Fractional 
errors in any quantities are less than 0.5%. 

As discussed in the section III C 5[ our choice of finite 
difference approximations is second order. The rate of 
convergence of our code is checked using different resolu- 
tions, whose setups of coordinate grids are shown in Table 
im The grid spacing of each coordinate (Ar, A9, A0) is 
proportionally scaled as 2/3, 3/4, 2/3, 3/4, from type 
I-l to T5. Here, we show the results of the convergence 
test with respect to the resolutions, fixing the maximum 



fl-i fc. 



AAl, + 0{Al+') 



(43) 



where fi-i (i = 1, • • • , 5) is a quantity computed using one 
of the resolution types 1-i in Table [Til /exact is its exact 
value, Ai-i represents the grid spacing associated with 
the type l-i setup, and ^ is a constant. Then, keeping 
the leading term, differences between different resolutions 
become 



/] 



i-fe 



Ai- 



- 1 



A^_ 



(44) 



To see the order n in a log-log plot, we select the com- 
binations of different resolutions that give the same ra- 
tio Ai-fc/Ai-j, and in our choice, these are /1-3 — /i-i, 
/1-4 - /1-2, and /1-5 - /1-3. In FigUl these combina- 
tions normalized by /1-5 of selected quantities are plot- 
ted against the grid spacing, where A represents the grid 
spacing in arbitrary units. It is clearly seen that the local 
quantities, here and e, converge to O(A^), and that 
integral (global) quantities also approach second order 
convergence as the resolution increases. 

We also checked the convergence with the different sets 
of resolutions type ITl, 2 and 3. These setups have fewer 
grid points in the angular coordinates 9 and (j). We found 
that the highest resolution type IT3 agrees well with the 
results of the higher resolutions of type 1-4 or 5 for the 
axisymmetric solutions. 



^ For the convergence tests with respect to the order of the Leg- 
endre expansion, see [Ssll 
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B. Triaxial solutions with M/R = 0.001 

We calculate triaxial sequences for small compactness 
M/R ~ 0.001, which is in the Newtonian regime, to 
check the value of at the bifurcation point of 

the triaxial sequence from the axisymmetric sequence. 
The triaxial and axisymmetric solution sequences for 
n = 0.3 and 0.5 are plotted in Fig. [2l Extrapolating 
the triaxial sequences to corresponding axisymmetric se- 
quences, values at the branch points are determined ap- 
proximately as [QM, T/\W\) = (1.735 x 10-^ 0.134) and 
(1.763 X 10"^ 0.135) for n = 0.3 and 0.5, respectively. 
This value may be compared with the Newtonian results 
such as the ellipsoidal approximation = 0.138 for 

n = 0.5 3. 



FIG. 1: The convergence of quantities, Q, A/adm, J, r/|W^| 
and e ~ \/l — (Rz/Rx)^ (in the proper length) for the 
model with n = 0.5, M/R = 0.296 and axis ratio in the 
coordinate length Rz/Rx = 0.75. Normalized differences 

l(/l-5 - /l-3)//l-5|, |(/l-4 - /l-2)//l-5|, and |(/l-3 " /l-l)//l-5| 

discussed in the text are plotted from left to right for each 
quantity against the resolutions A1-3, A1-2, and Ai-i, respec- 
tively. Black thin lines are proportional to A^. 
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FIG. 2: Plot of r/IWI versus normalized angular velocity 
nM for triaxial solution sequences (labeled by JB) are shown 
with curves marked with filled squares for n = 0.3, and with 
squares for n — 0.5. These sequences merge with axisym- 
metric solution sequences (labeled by ML) of the correspond- 
ing parameters (inset for a close up), which are shown by 
curves marked with plus (-|-) for n — 0.3, and with crosses 
(x) for n = 0.5. The compactness M/R of the sequences is 
set M/R = 0.001 for modeling the weak gravity regime. 



IV. TRIAXIAL SOLUTIONS 
A. Accuracy of tlie sequences of solutions 

Triaxially deformed solutions are calculated for se- 
lected values of the polytropic index, n — 0.3 and 0.5. 
Models with M/R ^ 0.001,0.1,0.14,0.2 are calculated 
for n = 0.3 and with M/R = 0.1,0.12,0.14 for n = 0.5. 

We noticed that it is necessary to increase the numbers 
of grid points as much as in type 1-5 in Table [TTl to have 
a smoothly changing sequence of triaxial solutions. For 
lower resolutions, the sequences appear to be less smooth 
especially for the plot of r/|W^| and for the part of the 
sequences closer to axisymmetric solutions. One of the 
reasons for this may be that, when one compares neigh- 
boring solutions of deformed sequences for mass, binding 
energy, or angular momentum, the change in these quan- 
tities for triaxial sequences is much smaller than that for 
axially symmetric sequences of about the same amount 
of deformation. 



B. Properties of the triaxial sequences 

As a sample of calculated solutions, the density con- 
tours in the xy, xz, and yz planes and the surface plot of 
the model with parameters n = 0.5 and M/R = 0.2 and 
the largest deformation are presented in Figl3l The solu- 
tion corresponds to the last row of data shown in Table 
El in Appendix [H 

In Figs. Hand El T/\W\ and nM are plotted for n = 
0.3 and 0.5 respectively, against the eccentricity for the 
constant rest mass sequences shown in the same Table 
IVl For uniformly rotating Newtonian polytropes, 
at the bifurcation point weakly depends on the difference 
of the EOS parameter, the polytropic index, whose value 
is about T/|iy| ~ 0.14. For highly differential rotations, 
T/\W\ may vary largely [H, [Hi]. The definitions of T 
and W in general relativity are given in Appendix [A] 

Our results for the solution sequences of uniformly ro- 
tating relativistic polytropes with n = 0.3 and 0.5 sug- 




DC 




x/R„ 



0.0 

y/Ro 





x/Rn 



FIG. 3: Contours of the p/p on a;z-plane (top left panel), 
on yz-plane (top right panel), and on a;?/- plane (bottom left 
panel) are shown for the most deformed triaxial model of n = 
0.3 and M/R = 0.2. Contours are drawn linearly from 0.0 to 
0.1 every 0.01 step. 



gest that the value of T/|H^| at the bifurcation point 
strongly depends on compactness M/R. In Tabic [TVl 
approximate values of quantities at the bifurcation point 
of each model are shown, which are evaluated by linearly 
extrapolating the triaxial sequence to the corresponding 
axisymmetric sequence. The value of T/|VF| at the bi- 
furcation point becomes ~ 0.169 for the compact model 
M/R = 0.2, n = 0.3, and it will certainly increase for a 
more compact sequence. 

As seen in the plot of Fig. [5l the triaxial solution se- 
quence for n = 0.5 becomes shorter as M/R increases. In 
fact, we were not able to find a triaxial solution sequence 
for M/R ~ 0.2; the triaxial sequence may disappear at 
a certain value of M/R between 0.14-0.2. We discuss an 
interesting consequence of the disappearance of triaxial 
sequences for high compactness in the last section. 



V. DISCUSSION: PROTO-NEUTRON STAR 
CONTRACTION 
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FIG. 4: Plots for T/\W\ (top panel) and QM (bottom panel) 
versus eccentricity e := \/l — {Rz/ RxY (iii proper length) 
for n = 0.3 sequences. Dashed curves labeled ML are ax- 
isymmetric solution sequences, and solid curves labeled JB 
triaxial solution sequences, where those correspond, from the 
top curves to the bottom in each panel, to M/R — 0.2, 0.14 
and 0.1 respectively. 



As a result of massive stellar core collapses, proto- 
neutron stars are formed and contract to more compact 
neutron stars within the time scale of cooling of a few tens 
of seconds [SB^ . Even for the small rotation rate of the col- 
lapsing stellar core, the ratio of the proto-neutron 
star becomes much higher than the value where the ax- 
isymmetric solution becomes secularly unstable against 
the viscosity driven £ = m = 2 bar mode instability 
[3, Q. Therefore, uniformly rotating triaxial solutions 
discussed in this paper may describe a quasi-stationary 
model of proto-neutron star contraction in the range of 
M/R ~ 0.1 — 0.2, assuming the following: (1) a cer- 
tain mechanism of strong viscosity operates during the 



contraction, (2) the time scale rapid cooling is shorter 
than that of gravitational radiation reaction, (3) the ef- 
fective polytropic (adiabatic) index n (F = 1 -I- 1/n) of the 
EOS for the realistic neutron star matter is small (large) 
enough to allow uniformly rotating triaxial solutions, and 
(4) those triaxial solutions are dynamically stable. 

Such an evolutionary track of a proto-neutron star con- 
traction has been considered using a compressible ellip- 
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n M/R R,/R:c £c » Madm J T/\W\ I Zp 

0.3 0.1 0.3226 (0.3718) 0.5634 (0.5693) 0.3888 0.6529 2.8444 x 10"^ 8.794 x 10"* 0.1507 1.347 x 10"^ 0.1328 

0.3 0.14 0.3435 (0.4220) 0.5531 (0.5619) 0.4445 0.7161 4.4203 x lO'^ 1.903 x 10"^ 0.1578 2.657 x 10"^ 0.2019 

0.3 0.2 0.3561 (0.4860) 0.5394 (0.5535) 0.5244 0.8063 6.9865 x 10"^ 4.351 x 10"^ 0.1688 5.395 x 10"^ 0.3311 

0.5 0.1 0.5153 (0.5919) 0.5461 (0.5536) 0.2059 0.4391 4.2103 x 10"^ 1.913 x 10"^ 0.1493 4.356 x 10"^ 0.1281 

0.5 0.12 0.5209 (0.6169) 0.5439 (0.5531) 0.2319 0.4698 5.2155 x 10"^ 2.753 x 10^^ 0.1519 5.861 x 10^^ 0.1594 

0.5 0.14 0.5314 (0.6482) 0.5366 (0.5482) 0.2581 0.4995 6.2304 x 10"^ 3.734 x 10'^ 0.1542 7.476 x 10'^ 0.1932 

TABLE IV: Quantities at the point of bifurcation of triaxial sequences from axisymmetric sequences. The polytropic index n 
and the compactness of the spherical star with the same rest mass M/R are the model parameters. Corresponding triaxial 
sequences are found in Table |V]in Appendix[B] In the above, Rx is the equatorial radius, and Rz/Rx is the ratio of polar to 
the equatorial radius. Each has two values; one is measured in the coordinate length, and the other in parenthesis is in proper 
length defined in Ea. ()AlH ). tc is the energy density at the center of the compact star, is the angular velocity. Definitions of 
Madm, J , T /\W\, and / are found in Appendix|X] Zp is the polar redshift. Dimensional quantities are shown inG = c = fi: = l 
units. 
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soidal model^ [H] ■ Our results add two further important 
features to this. First, the sequences of triaxial solu- 
tions terminate at the maximally deformed models, at 
the mass-shedding limits, and the changes in ADM mass 
or total angular momentum are small along the triaxial 
sequences from the bifurcation points to the termination 
points, even for the stifFer EOS such as n = 0.3, as seen 
in Table |Vl the triaxial sequences are not very long at 
all. Secondly, the triaxial sequences may become shorter 
and disappear as the compactness becomes larger for a 
relatively less stiff EOS such as n = 0.5. 

The angular velocity near the braking limit is esti- 
mated as nM ~ {M/Rf^ - {M/Rfj/M^. Therefore 
once the secular bar mode instability sets in, conserving 
M and J, the proto-neutron star evolves towards a maxi- 
mally deformed triaxial configuration as it contracts, say, 
from M/R ~ 0.1 to 0.2. And then, it is likely that it al- 
ways evolves along the sequence of maximally deformed 
configurations during the contraction as long as such tri- 
axial solutions exist and are dynamically stable in the 
parameter region of the effective T and M/R. Excess an- 
gular momentum arising from contraction may be trans- 
ported outward by mass ejected from the Lagrange point 
at the cusp of the longest semi-major axis. (Note that 
the time scale of the mass ejection may be that of cool- 
ing, which is much longer than the dynamical time scale.) 
Furthermore, if the effective F satisfies F < 3, the triaxial 
solution may disappear when the solution reaches a cer- 
tain value of the compactness M/R and higher. Dynam- 
ically stability of such uniformly rotating solutions are 
not known, but it is unlikely that the dynamical insta- 
bility appears within such short triaxial sequences, along 
which the ratio is nearly constant. 

It is estimated that the amplitude of gravitational wave 
(GW) signals from such objects may be detectable using 
the ground based laser interferometric detectors, if the 



FIG. 5: Same as Fig. |4] but for n = 0.5 sequences. Dashed 
curves and solid curves from the top to the bottom in each 
panel correspond to M/R = 0.14, 0.12 and 0.1 respectively. 



In their work, changes in the entropy during the evolution is mod- 
eled by the changes in the adiabatic constant k of the polytropic 
EOS. 
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source is within a few tens of Mpc [1, [T^ . Detection of the 
persistent GW signals even after the proto-neutron star 
contraction phase suggests a large effective F > 3, while 
the shutdown of the signal during the contraction implies 
the relatively smaller F < 3. Detection of such GW sig- 
nal may set another constraint on the EOS parameter of 
high density matter. Source modeling for constructing 
the wave templates may be straightforward because one 
can concentrate on calculating the maximally deformed 
configurations. Our next plan is to include more realis- 
tic nuclear EOS in the code, then to estimate the grav- 
itational wave amplitude for those EOSs that allow the 
triaxial solutions. 
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APPENDIX A: FORMULAS FOR MASS AND 
ANGULAR MOMENTUM 

Definitions of the quantities shown in tables and figures 
that characterize each solution of a rotating relativistic 
star, and their expressions in terms of the metric poten- 
tials in the IWM formulation, are summarized in this 
Appendix. 

The rest mass of the star Mq is written as 



Mq / pu"dSa = / pu^a^P^y/f(fx (Al) 

where dSa — at\/—gcPx and ^/—gdPx — aip^y^cfix — 
at/j^r'^ sin ddrd6d(f>. 

The ADM mass Madm becomes 



Ma 



DM 



1 f (f-^fd-f-bfcd^^D^^^^dSa 

167r 



2n 



D'^^dSa 



2n 



Ail^dS 



2-K 



1 



(A2) 



where dSa — ar\fjd'^x and dS = y/fd'^x, and dSa 

o 

coincides with dSa at spatial infinity. 

The Komar mass associated with a timelike Killing 
field t" is written 

Mk - -J^J ^"f'dSc^p = -1. J^R'-f^t" dSa 



j^{2T''p-Tgy)t^ dS^, 



(A3) 



and, in the IWM formulation, we have 

Mk= / [a{pii + S)-2jaf3'']^^y^d^x, (A4) 

where dSa — n.ay^d'^x was used. The above derivation 
holds if the global timelike Killing field exists. For the 
spacetime of a triaxially deformed rotating star, no such 
timelike Killing field exists. Instead, an asymptotic Ko- 
mar mass can be written 



Mk := - 



1 



47r 



47r 



= — AadT, 



aiabi"'' + 4^a (pH + S)j ^p'^y/fd^x. 

(A5) 



In [35|, we have derived sufficient conditions of the fall 
off of the 3-metric 7ab and extrinsic curvature Kat and 
their time derivative for the Madm = Mk relation to be 
satisfied. In the IWM formulation the fall off of each field 
is sufficiently fast to have the equality. And also in this 
case, the above two definitions for Mk agree. 

The total angular momentum calculated in the asymp- 
totics is written 



60") dS 



(A6) 
(AT) 



The relativistic analog of the kinetic energy T is de- 
fined by 



(A8) 



therefore, for uniform rotation we have T = \nj . Also 
the relativistic analog of the gravitational potential en- 
ergy W is defined by 



I^:=A/p + r-MADM, 
where Mp is the proper mass defined by 



Mn 



eu'^dSa 



eu^a 



i^^^/fd'x. 



(A9) 



(AlO) 



The proper length of the semi-major axis along the x 
direction is written 



R:r. ^ I 1p dx, 



(All) 
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where is the coordinate length of the same axis. The 
proper lengths along the y or z directions are calculated 
using the same formula, replacing xhy y oi z respectively. 

In the above, the source terms of the field equations, 
Ph, Ja, and S*, are obtained from the stress energy ten- 
sor. We write down the projection of the stress tensor 
in terms of elementary fluid variables and metric poten- 
tials. The 4-velocity for the corotational flow u" = m'/c" 
is decomposed with respect to the foliation Et as 

M"na = -au* (A12) 
u"7qq = u'wa- (A13) 

Using these relations, the source terms of the field equa- 
tions become 

Ph Ta,0n°'n^ = hp{au* f ~ p, (A14) 
Ja -Taf^la^'n^ = hpa{u'f^^Cja. (A15) 
S := T^p-f'^^ ^ hp[{au' f -l\+?,p, (A16) 

where Ua /ab^' = fabiP'' + f^^') = Pa + ^k- 
Throughout the paper we use units such that G ~ c ^ 



K = 1. The latter equality is implemented by renormal- 
izing the length and mass scales as 

R M := k.-'^/'^M. (A17) 

respectively. Angular momentum and angular frequency 
are respectively normalized as 

J K-"J, n := (A18) 

We omit the bars over these quantities in the main text. 



APPENDIX B: SELECTED SOLUTION 
SEQUENCES 

Sequences of triaxially deformed solutions of the com- 
pact star were calculated for the following parameters: 
the polytropic index n — 0.3 with compactness M/R = 
0.1, 0.14, 0.2, and the polytropic index n = 0.5 with com- 
pactness M/R — 0.1,0.12,0.14. Quantities, most of 
which are defined in Appendix [X] are tabulated in Ta- 
ble El This data is plotted in Figlland[5] 
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n = 0.3 


M/i? = 0.10 


Mo = 


: 2.9908 X 10"^ M = 


-- 2.8116 X 10" 


■2 








Rx 


Ry/ Rx 


Rz/ Rx 






Madm 


J 




T/\W\ 


I 




Zp 


0.3367 (0.3875) 


0.9375 (0.9389) 


0.5423 (0.5485) 


0.3884 


0.6504 


2.8446 X 10"^ 


8.845 X 10' 


-4 


0.1512 


1.360 X 10" 


-3 


0.1327 


0.3461 (0.3981) 


0.8906 (0.8929) 


0.5273 (0.5337) 


0.3882 


0.6493 


2.8447 X 10"^ 


8.871 X 10- 


-4 


0.1515 


1.366 X 10" 


-3 


0.1326 


0.3567 (0.4098) 


0.8438 (0.8470) 


0.5115 (0.5181) 


0.3880 


0.6475 


2.8449 X 10"^ 


8.909 X 10- 


-4 


0.1518 


1.376 X 10" 


-3 


0.1326 


0.3645 (0.4185) 


0.8125 (0.8164) 


0.5008 (0.5075) 


0.3878 


0.6457 


2.8450 X 10"^ 


8.937 X 10" 


-4 


0.1520 


1.384 X 10" 


'3 


0.1325 


0.3775 (0.4330) 


0.7656 (0.7704) 


0.4833 (0.4902) 


0.3874 


0.6423 


2.8453 X 10"^ 


8.992 X 10' 


-4 


0.1524 


1.400 X 10" 


^3 


0.1323 


0.3875 (0.4441) 


0.7344 (0.7397) 


0.4707 (0.4778) 


0.3871 


0.6397 


2.8455 X 10"^ 


9.043 X 10- 


-4 


0.1528 


1.414 X 10" 


-3 


0.1322 


0.3987 (0.4564) 


0.7031 (0.7091) 


0.4574 (0.4647) 


0.3867 


0.6366 


2.8458 X 10"^ 


9.095 X 10" 


-4 


0.1532 


1.429 X 10" 


-3 


0.1320 


0.4188 (0.4784) 


0.6562 (0.6633) 


0.4355 (0.4433) 


0.3863 


0.6318 


2.8461 X 10"^ 


9.171 X 10' 


-4 


0.1536 


1.452 X 10" 


-3 


0.1318 




n = 0.3 


M/R = 0.14 


Mo = 


: 4.7471 X 10"^ M = 


4.3417 X 10" 


-2 








Rx 




Rz/ Rx 






Madm 


J 




T/\W\ 


/ 







0.3619 (0.4435) 0.9219 (0.9245) 0.5291 (0.5385) 0.4441 0.7124 4.4208 x 10"^ 1.913 x 10"^ 0.1581 2.686 x 10"^ 0.2017 

0.3763 (0.4603) 0.8594 (0.8639) 0.5089 (0.5187) 0.4439 0.7103 4.4210 x 10"^ 1.919 x 10"^ 0.1582 2.702 x 10~^ 0.2016 

0.3889 (0.4751) 0.8125 (0.8185) 0.4924 (0.5026) 0.4436 0.7077 4.4213 x 10"^ 1.927 x 10"^ 0.1584 2.722 x 10"^ 0.2014 

0.3985 (0.4861) 0.7812 (0.7882) 0.4807 (0.4911) 0.4434 0.7056 4.4216 x 10"^ 1.933 x 10"^ 0.1586 2.739 x 10"^ 0.2012 

0.4091 (0.4984) 0.7500 (0.7579) 0.4684 (0.4792) 0.4431 0.7032 4.4220 x 10"^ 1.940 x 10"^ 0.1587 2.758 x 10"^ 0.2010 

0.4215 (0.5127) 0.7188 (0.7278) 0.4548 (0.4659) 0.4429 0.7007 4.4223 x 10"^ 1.947 x 10"^ 0.1589 2.778 x 10"^ 0.2008 

0.4364 (0.5296) 0.6875 (0.6978) 0.4395 (0.4513) 0.4427 0.6983 4.4226 x 10'^ 1.952 x 10"^ 0.1590 2.796 x 10'^ 0.2007 

n = 0.3 M/R = 0.20 Mq = 7.7530 x 10'^ M = 6.7804 x 10"^ 

Rx Ry/Rx Rz/Rx £c n A/adm J T/\W\ I Zp 

0.3745 (0.5091) 0.9219 (0.9262) 0.5162 (0.5313) 0.5242 0.8036 6.9868 x 10"^ 4.360 x 10"^ 0.1687 5.426 x 10"^ 0.3308 
0.3819 (0.5184) 0.8906 (0.8967) 0.5064 (0.5219) 0.5242 0.8026 6.9869 x 10"^ 4.362 x 10"^ 0.1686 5.435 x 10"^ 0.3306 
0.3945 (0.5343) 0.8438 (0.8524) 0.4905 (0.5066) 0.5241 0.8008 6.9870 x 10"^ 4.369 x 10"^ 0.1686 5.456 x 10"^ 0.3305 
0.4227 (0.5691) 0.7656 (0.7791) 0.4584 (0.4762) 0.5240 0.7968 6.9873 x 10"^ 4.380 x 10"^ 0.1683 5.498 x 10'^ 0.3300 

n = 0.5 M/i^ = 0.10 Mq = 4.4113 X 10'^ M = 4.1580 x 10'^ 

Rx Ry/Rx Rz/Rx ec n Madm J T/\W\ I Zp 

0.5407 (0.6197) 0.9219 (0.9243) 0.5238 (0.5319) 0.2059 0.4380 4.2103 x 10"^ 1.915 x 10"^ 0.1492 4.373 x 10"^ 0.1281 
0.5580 (0.6387) 0.8750 (0.8789) 0.5079 (0.5163) 0.2059 0.4373 4.2104 x 10"^ 1.917 x 10"^ 0.1491 4.384 x lO"'' 0.1280 
0.5716 (0.6536) 0.8438 (0.8486) 0.4959 (0.5047) 0.2059 0.4367 4.2104 x 10"^ 1.919 x 10"^ 0.1491 4.393 x 10"^ 0.1280 
0.5880 (0.6715) 0.8125 (0.8185) 0.4823 (0.4914) 0.2059 0.4361 4.2104 x 10"^ 1.920 x 10"^ 0.1490 4.402 x lO"'^ 0.1279 

n = 0.5 M/R = 0.12 Mp = 5.5171 x 10'^ M = 5.1345 x 10'^ 

Rx Ry/Rx Rz/Rx £c » A/adm J T/\W\ I Zp 

0.5508 (0.6500) 0.9219 (0.9250) 0.5188 (0.5289) 0.2319 0.4689 5.2155 x 10"^ 2.755 x 10"^ 0.1517 5.875 x 10"^ 0.1594 
0.5693 (0.6708) 0.8750 (0.8800) 0.5021 (0.5127) 0.2319 0.4683 5.2156 x 10"^ 2.758 x 10"^ 0.1517 5.888 x 10"^ 0.1594 
0.5934 (0.6976) 0.8281 (0.8352) 0.4818 (0.4932) 0.2319 0.4677 5.2156 x 10"^ 2.759 x 10~^ 0.1516 5.900 x 10"^ 0.1593 

n = 0.5 M/R = 0.14 Afp = 6.6547 x 10"^ M = 6.1130 x 10"^ 

Rx Ry/Rx Rz/Rx £c n A/adm J T/\W\ I Zp 

0.5410 (0.6591) 0.9688 (0.9703) 0.5282 (0.5401) 0.2580 0.4993 6.2304 x 10"^ 3.735 x 10"^ 0.1542 7.482 x 10"^ 0.1932 
0.5630 (0.6843) 0.9062 (0.9109) 0.5077 (0.5203) 0.2580 0.4988 6.2305 x 10"^ 3.738 x 10"^ 0.1542 7.495 x 10"^ 0.1931 
0.5771 (0.7003) 0.8750 (0.8813) 0.4954 (0.5085) 0.2580 0.4984 6.2306 x 10"^ 3.740 x 10"^ 0.1542 7.505 x 10"^ 0.1931 

TABLE V: Quantities for the constant rest mass Jacobi-like triaxial sequences. Mp is the rest mass of each sequence, and A/ is 
the gravitational mass of the spherical star having the same Mo. In the definition of the compactness, R is the Schwartzschild 
radius of the spherical star. Tabulated quantities are the same as Table ITVl except for Rx, which is the semi-major radius along 
the X-axis, and Ry, the semi-major radius along the a;-axis. Dimensional quantities are shown inG = c = K= l units. 



